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Abstract 

A model describing spin-polarized current via discrete energy levels of a metallic nanoparticle, 
which has strongly asymmetric tunnel contacts to two ferromagnetic leads, is presented. In absence 
of spin-relaxation, the model leads to a spin- accumulation in the nanoparticle, a difference (A//) 
between the chemical potentials of spin-up and spin-down electrons, proportional to the current and 
the Julliere's tunnel magnetoresistance. Taking into account an energy dependent spin-relaxation 
rate rt[ijo), Afi as a function of bias voltage (V) exhibits a crossover from linear to a much weaker 
dependence, when |e|ri(A/i) equals the spin-polarized current through the nanoparticle. Assuming 
that the spin-relaxation takes place via electron-phonon emission and Elliot- Yafet mechanism, 
the model leads to a crossover from linear to V^^^ dependence. The crossover explains recent 
measurements of the saturation of the spin-polarized current with V in Aluminum nanoparticles, 
and leads to the spin-relaxation rate of ~ 1.6MHz in an Aluminum nanoparticle of diameter 6nm, 
for a transition with an energy difference of one level spacing. 



INTRODUCTION 



Spin-dependent electron transport through nanometer scale structures has attracted an 
increased interest recently, [l, 2] In general, short spin-diffusion lengths in metals makes it 
necessary to investigate spin-dependent transport in micron-scale metallic structures. [3|, |4| 
More recently, spin-dependent electron transport has been investigated in single nanometer- 
scale metallic particles, 5j_^ including spin-polarized transport via discrete electronic energy 
levels of a nanoparticle. 6| 

In this paper, we develop a model to explain a previously reported experiment designed 
to detect spin polarized currents in a normal Aluminum nanoparticle connected to ferromag- 
netic leads by weak tunnel barriers. These experiments were carried out on lithographi- 
cally defined tunnel junctions as featured in fig. [H In principle, spin polarized current can 
be determined by analyzing the TMR = (/|| — I-^i)/!-^-^, where J|| and are the currents 
through the nanoparticle in the parallel and antiparallel magnetization configurations. 

It is assumed here that the difference between and I^i arises from the differences 
in spin-up and spin-down densities of states in the leads. In sequential electron tun- 
neling via the nanoparticle, spin-polarized current is a consequence of spin accumula- 
tion, a difference in the chemical potentials of the spin-up and spin-down electrons in the 

electric current and the spin-polarized densities of states. 
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in the antiparallel magnetization configuration, because only in that case the ratios of the 
tunnel- in and tunnel-out resistances of the spin- up and the spin down electrons are different. 
A necessary condition for spin-accumulation is that the electron spin be conserved during 
the sequential transport process. 

There are several compelling reasons to study TMR in nanometer scale particles at tem- 
peratures where discrete energy levels can be resolved. One is the magneto-coulomb effect 
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231] which can give a strong TMR signal even without spin-accumulation in the nanopar- 



ticle. In the regime of well resolved energy levels the contribution arising from spin-polarized 
current and spin-accumulation can be separated from the contribution arising from the chem- 
ical potential shifts. 6| In particular, if the leads chemical potentials vary in a range that 
corresponds to a nanoparticle energy range that is in between two successive discrete energy 
levels, then electron current through the nanoparticle will be insensitive to changes in the 
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chemical potentials. In terms of the I-V curve, which exhibits step like increases at bias volt- 
ages corresponding to discrete energy levels, TMR must be measured between the current 
steps, at voltages where the current as a function of bias voltage is constant. 

Another reason is that nanoparticles of this size exhibit extraordinarily weak spin-orbit 



coupling compared to bul 
perturbed spinors. 
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i. In this regime, the stationary electron wavefunctions are slightly 



25l | As a result, the spins of electrons injected from a ferromagnet into 



the nanoparticle have exceptionally long life-times. This remarkable regime of spin-polarized 
electron transport via metallic nanoparticles has hardly been explored experimentally. 5, 6| 
By contrast, most measurements of electron spin-injection and accumulation in metals 3|, |4| 
were obtained in the regime of strong spin-orbit scattering, where any nonequilibrium spin- 
population exhibits exponential time-decay. 

The main result of the prior report {ff] was that the spin-polarized current through the 
nanoparticle saturates quickly as a function of bias voltage, for tunnel resistances in the Gfi 
range. In that case the saturation is reached typically around the second or third energy 
level of the nanoparticle. 

The saturation effect was explained by a rapid increase of the spin-relaxation time with 
the nanoparticle excitation energy. By this interpretation, the spin-polarized current through 
the nanoparticle is carried only via the ground state and the few lowest energy excited states 
of the nanoparticle, while highly excited spin-polarized states relax faster than the average 
sequential electron tunnel process. We conjectured that at the saturation voltage, the re- 
laxation time of the highest singly occupied energy level of the nanoparticle is comparable 
to the electron tunnel rate. The corresponding spin relaxation time is in the microsecond 
range. In comparison, the spin relaxation time in Aluminum thin films with a similar mean 
free path would be five orders of magnitude shorter. 

Our measurement of the spin-relaxation time was somewhat indirect, because theoretical 
literature prior to our work had not predicted any saturation of the spin-polarized current 



with bias voltage. 
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221] Consequently 



our explanation of the saturation was qualitative. The goal of this paper is to obtain a 
model of spin-polarized electron transport through a metallic nanoparticle to explain our 
observations. We use a method for calculating current via energy levels of the nanoparticle 
based on rate equations, following Ref . [26i] . The tunneling regime in our devices is different 
from that used in the theoretical studies, because the tunnel resistances in our junctions are 
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FIG. 1: A: Sketch of sample cross-section in the antiparahel magnetization configuration of the 
leads. B. Scanning Electron micrograph of a typical device. The tunneling junction with embedded 
nanoparticles is located inside the overlap of two Py leads. 



highly asymmetric and spin-relaxation rate has strong energy dependence 
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Other experimental work on arrays and single nanoparticles did not find any saturation 
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1111] These experiments do not 



of the spin-polarized current with bias voltage. [5|, 
measure TMR in the regime of well resolved energy levels, we believe that this is a critical 
measurement to separate the contributions to TMR from the chemical potential shifts. The 
analysis of TMR experiments in Ref. ||1|] uses an energy independent spin-relaxation time. 
It is possible, however; that the energy dependence of the spin-relaxation rate can be quite 
strong. 291] We show in this paper that the effect of the energy dependence on spin-polarized 
current is significant. 

Our model is valid within a specific experimental regime, outlined in Sec. , but it is an 
easily analyzable and experimentally relevant regime that merits consideration as a large 
number of samples fall under this parameter range. 

In Sec. we review the effects of spin-orbit scattering in metallic nanoparticles, in the 
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context of spin-injection and detection. In Sec. we discuss various energy-relaxation rates 
in the nanoparticle. The criteria of the model validity are listed on Sec. . We then calculate 
the probability distribution of the many-electron states in the nanoparticle in Sec. and use 
this to calculate a TMR versus bias voltage curve that can be fit to experimental data in 
Sec. 



EFFECTS OF SPIN-ORBIT INTERACTION ON SPIN-POLARIZED CURRENT 
THROUGH A NANOPARTICLE. 

In a metallic nanoparticle, the stationary electronic wavefunctions in zero applied mag- 
netic field are two fold degenerate and form Kramers doublets: 

I ]>'=u{r)\ \>+v{r)\ i>, 

||>'=n^(f)|i>-.^(f)|T>. 

The mixing between the spin-up and the spin-down components is caused by the spin-orbit 
interaction. In a magnetic field, the degeneracy is lifted by the Zeeman effect. If the field is 
applied in direction corresponding to | ]'>, the g-factor is = 2 — 4 / (i^r|f p < 2. 

The effects of the spin-orbit interactions on the wavefunctions are described by a dimen- 
sionless parameter A 

A = M^, (1) 



where vso is the spin-orbit scatterin 



Kramers doublets in the particle. 



I rate and 5 is the average spacing between successive 
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30| The effects of spin orbit scattering are weak if 
A <^ 1. In that case, v{r) ^ 0, u{f) is real, and the wavefunctions | t>' and | |>' have well 
defined spin. In the opposite limit, A ^ 1, the spin-orbit scattering is strong. In that case, 
/ (i^rl-up ^ J (i^r|f p 1/2 and the wavefunctions have uncertain spin. 

Consider an electron with spin-up injected by tunneling from a ferromagnet into the 
nanoparticle. If the tunnel process is instantaneous, the electron will have a well defined 
spin I ]'> immediately after tunneling. In the regime of weak spin-orbit scattering, where 
f ~ 0, the initial state has a much larger overlap with a wavefunction | ]■>' than with the 
corresponding wavefunction | |>'. As a result, the spin of the added electron will remain 
well defined, in principle indefinitely long, barring any coupling between the nanoparticle 
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and the environment. In that case, the detection of the injected spin can be performed at 
any time after the injection. 

By contrast, in the regime of strong spin-orbit scattering, the initial state overlaps with 
both I ]■>' and | ],>', nearly equally. In that regime, the spin of the added electron becomes 
uncertain after a time ~ tso = ^/^SO) so there is a time limit ~ tso for spin detection. 

In normal metals, there is a scaling between the spin-conserving electron scattering rate 



and the corresponding spin-flip electron scattering rate. 



^so{'^) = au{uj). 
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(2) 



This is known as the Elliot- Yafet relation, u is the energy difference between the initial 
and final state. For elastic scattering, u = and 1^50(0) = i^so- The scaling parameter 
a depends on the spin-orbit scattering and the band structure. In Aluminum, a ~ 10~^ 
is larger than anticipated from the spin-orbit interaction, because of the hot-spots for spin 
scattering in the band structure. |32l. l33| 

The Elliot-Yafet relation was confirmed in bulk metals by the conduction-electron spin 
resonance experiments (CESR). In particular, the temperature dependence of the width 
of the spin-resonance line, which is proportional to usoiknT), follows the temperature 
dependence of the resistivity, which is proportional to the momentum relaxation rate, in 



agreement with Eq. [2J The relation has also been confirmed more recent 
metallic samples, by the spin-injection and detection experiments. {4, 34, 



y m mesoscopic 



35j Both CESR 



and spin-injection and detection experiments measure the time decay of a nonequilibrium 
spin population in the regime of strong spin-orbit scattering. 

In metallic nanoparticles, the validity of the Elliot-Yafet relation was confirmed exper- 
imentally by energy level spectroscopy, in the regimes of weak and moderate spin-orbit 
scattering. 13611 uso can be measured directly from the magnetic field dependence of the en- 
30 1 Petta et al. found that in Cu, Ag, and Au nanoparticles, vso 



ergy levels. 



where vp is the Fermi velocity and D is the nanoparticle diameter, and a is close to the bulk 
value, within an order of magnitude. This confirms the Elliot-Yafet relation because vpjD 
is equal to the elastic scattering rate, assuming a ballistic nanoparticle. More recently, the 
Elliot-Yafet relation was confirmed in Al nanoparticles as well. 6| 

Substituting the level spacing and the Elliot-Yafet relation into Eq. [T|, we find that there 
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is a characteristic nanoparticle size 



= Xf/V^. (3) 

Spin-orbit scattering will be weak if D ^ D* and spin-orbit scattering will be strong if 

45l | In Aluminum, ^ 



D ^ D* . D* is a material dependent microscopic parameter. 
IO^Af ~ lOnm. 

It follows that both mesoscopic and macroscopic metallic samples are in the regime of 
strong spin-orbit scattering. Only if the nanoparticle diameter is less than about lOnm 
the spin-orbit scattering becomes weak. In samples much larger than D*, measuring the 
discrete electron energy levels and the electron spin polarization are incompatible, and the 
spin-polarized current via resolved energy levels must be negligibly small. We are not aware 
of any theory that explicitly takes into account the effect of spin-mixing in the Kramers 
doublets on TMR of ferromagnetic single-electron transistors, e.g. that calculates TMR 
versus A. Such a theory would be important because it would set the limits of observability 
of spin-polarized current through single electron transistors. 

In samples much larger than D*, one can still study the time decay of injected electron 
spins at a time scale much shorter than the time necessary to resolve discrete energy levels 
(the Heisenberg time h/6), as demonstrated by the CESR experiments and spin-injection 
and detection in mesoscopic metals. In further discussion, we assume D ^ D*. In that case, 
the spin up band in the ferromagnet is tunnel coupled to the nanoparticle wavefunctions 
with spin-up only. 



ENERGY RELAXATION IN METALLIC NANOPARTICLES 

In this work we study spin-polarized electron current via discrete energy levels of metallic 
nanoparticles, in a regime where the tunnel rate is much smaller than the spin-conserving 
energy relaxation rate. The tunnel rate can still be larger than the spin-flip relaxation rate, 
so electrons in the nanoparticle are not in equilibrium. 

Consider a nanoparticle with an electron added at energy u = 6, a.s shown in Fig. [2l The 
nanoparticle is in an excited state and it can relax its energy internally. It was shown that 
the dominant relaxation process is electron-phonon interaction, js?] One possible path for 
this relaxation is spin-conserving, as shown in Fig. [2]-A. The spin-conserving relaxation rate 
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FIG. 2: A: Spin-conserving energy relaxation process. B. Spin-flip energy relaxation process. C: 
Electronic state with the relaxation time equal to the spin- flip relaxation time. 



via phonon emission have been estimated by Agam et al. [3 



(4) 



where Ep = 11.7eV is the Fermi energy, u is the energy difference between the initial and 
the final state, is the elastic scattering relaxation time, p = 2.7g/cm^ is the ion-mass 
density, and Vs = 6420m/ s is the sound velocity. This formula is equivalent to 



hVe-ph{'^) 



g \^dJ 



within a prefactor of order 1, where tu^ is the Debye energy and g is the dimensionless 
conductance of the nanoparticle, g = Eth/S, and Eph = hvp/D is the Thouless energy. This 
relation is similar to the formula for bulk electron-phonon scattering rate, except that there is 
a prefactor 1/g, [37| which originates from the chaotic nature of the electron wavef unctions. 
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We estimate Ve-ph{8) ~ GHz in nanoparticle with diameter 5nm, in agreement with 



experiment. 



38l | All our samples have tunnel rates significantly smaller than this relaxation 



rate, so, the nanoparticle can be considered to be relaxed toward the lowest energy state 
accessible via spin-conserving transitions. 

The nanoparticle can also relax through a transition depicted in Fig. [2]-B, with some 
rate v^^{uj)^ uo = 5. This spin-flip transition involves coupling between electrons and the 
environment, which may be the phonon bath or the bath of nuclear spins. In zero applied 
magnetic fleld, the electron transitions shown in Figs. [2]-A and B have the same energy 
difference between the initial and the flnal states. We expect that the transition rate for the 
spin-flip relaxation process is much smaller than the corresponding spin-conserving relax- 
ation rate, because in the spin-flip process there must be a transfer of angular momentum 
into the environment. 

For example, if the environment is the ion lattice, then the transfer is governed by the 
spin-orbit interaction, which is much smaller than the the electrostatic interaction that 
governs the spin-conserving transitions. Thus the spin-flip probability caused by the spin- 
orbit interaction during a phonon emission process is very small. It is reasonable to expect 
that the Elliot- Yafet scaling is valid for the transitions between chaotic wavef unctions, so 
that v^^{uj) = ave-phi'^)-, but we are not aware of any theoretical calculation of the spin 
relaxation rate in metallic nanoparticles. 

Spin-flip electron transition rates between discrete levels were obtained theoretically for 
semiconducting quantum dots. 39|, |40[ The theoretical calculations were in good agreement 
with the experimental results in GaAs quantum dots. 4l| The theory of spin-relaxation in 
semiconducting quantum dots is of no use for metallic nanoparticles, because the mechanisms 
of spin-orbit interaction in metals and semiconductors are very different. [? ] 

It would be very difficult to measure the transition rate depicted in Fig. [2]-B, because the 
relaxation process in Fig. [2]-A competes with that in Fig. [2]-B. We are able to determine the 
spin relaxation rates in tunneling measurements because it is possible to trap the nanopar- 
ticle into a state shown in Fig. [2]-C. In that case the spin-conserving energy relaxation is 
forbidden by the Pauli principle and the relaxation rate is equal to the spin relaxation rate. 
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REGION OF MODEL VALIDITY 



In this paper the tunnel current via discrete energy levels of a metallic nanoparticle is 



calculated using the methods outlined in Ref. [26j We consider spin polarized tunneling in 
the regime closest to our experiments. The following conditions define that regime. 

1. We assume that the tunnel junctions in our samples are highly asymmetric in resis- 
tance. The asymmetry arises because the tunnel junctions are fabricated using conventional 
lithography and evaporation. A weak non-uniformity in the tunnel junction thickness leads 



to large asymmetry in the tunnel resistance. 



6| 



In these asymmetric samples, the current through the nanoparticle is limited by the tunnel 
rate through the high resistance junction. The resistance of the low resistance junctions is 
less directly related to the current and it can be obtained by comparing the amplitudes of 
Zeeman split energy levels at positive and negative bias voltage. That procedure has been 



described in detail in some special cases in Ref. [26| Even though our regime is quite different 
from those special cases, they are still indicative of the procedure to determine the resistance 
ratio. Applying the general formalism of Ref. to our regime we estimate the resistance 
ratio to be about 25 in sample 1. The results of our model do not depend on the resistance 
ratio, as long as the electron discharge rate is much larger than the electron tunnel in rate. 

2. We calculate the current in the regime when the first tunnel step is across the high 
resistance junction. This step is followed by an electron discharge via the low-resistance 
junction. This regime is relatively easy to analyze because electron discharge is fast and 



so the charging effects do not influence the current significantly. [26|, |42] In this case, the 
nanoparticle spends most of the time waiting for an electron to tunnel in. 

3. We assume Ec ^ S, where Ec is the charging energy. This assumption is generally 
valid in metallic nanoparticles. 38 1 

4. We assume that the coulomb gap in the I-V curve is much larger than the level 
spacing, this requires not only that condition 3 is met, but that the background charge go 
is not too close to {n + l/2)e, where n is an integer. In that case the number of energy 
levels participating in electron transport is always ^ 1. Even if an electron tunnels into the 
lowest unoccupied single-electron energy level of the nanoparticle, there will still be a large 
number of occupied single-electron energy levels of the nanoparticle that can discharge an 
electron. 
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5. We assume that the number of electrons on the nanoparticle before tunnehng in is 
even. The calculation for the odd case is very similar to that for the even case and will not 
be discussed here. 

6. We assume that the tunnel rates between the single-electron states i in the nanoparti- 
cle and the leads are much smaller than the spin- conserving energy relaxation rates, which 
are > GHz. In the sample selected for this paper, the tunnel rate across the high resistance 
junction is in the MHz range. 

If the spin-relaxation process is taking place and there is a large asymmetry in junction 
resistances, an important question is which one of the two tunnel rates limits the TMR. We 
will assume that the left junction has higher resistance. It may be tempting to assume that 
TMR will be highly asymmetric with bias voltage if -C u^^luj) <^ F^, where F^ and 
F^ are the electron tunnel rates between the energy levels of the nanoparticle and the left 
and the right leads, respectively, because, if a spin-polarized electron tunnels in via the high 
resistance junction, it will tunnel out via the low resistance junction before spin-relaxation 
process takes place, and, if the direction of the current is reversed, the order of tunneling will 
be reversed and the spin-relaxation will take place before tunneling out. A similar situation 
is found in measurements of the energy spectra in samples where F^ <^ u^-phi^) ^ T^j 
where i>e-ph{uj) is the spin-conserving energy relaxation rate. In that case, electron trans- 
port is much closer to equilibrium in one direction of current than in another direction of 
the current, resulting in asymmetric energy level spectra, jssi 




It turns out that in spin injection and detection in the regime defined in this section, 
TMR is symmetric even if F^ <C u^^^u) <^ F^. The reason is that the spin-polarized cur- 
rent is mediated by spin-accumulation, which takes place after a large number of tunnel-in 
and tunnel-out steps. 

Assume that an electron first tunnels in via the low resistance junction and then an elec- 
tron tunnels out via the high resistance junction. In that case, it is clear that in order to 
observe spin polarized current, it is necessary that the spin-relaxation time in the nanopar- 
ticle be longer than the tunnel out time. 

If the bias voltage is reversed, an electron first tunnels in via the high resistance junction 
and then an electron tunnels out via the low resistance junction. In our regime, it is highly 
improbable that the same electron tunnels in and out, because the number of occupied elec- 
tron states available for discharge is ^ 1. To obtain a spin-polarized current in this case, the 
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spin-accumulation in the nanoparticle is necessary, which takes place after many tunnel in 
and tunnel out steps and makes it necessary that the spin of the nanoparticle be conserved 
during the time that the nanoparticle waits before an electron tunnels in. Overall, the spin- 
polarized current and TMR are comparable in magnitude for the two current directions, in 
agreement with our measurements. 



CALCULATION OF THE SPIN PROBABILITY DISTRIBUTION IN THE 
NANOPARTICLE 



In this section we obtain the probability distribution among various many-electron states 
I a > generated by electron tunneling via the nanoparticle. The many-electron states are 
the Slater determinants of varying single-electron states of the nanoparticle. In further 
discussion, the many electron states will be referred to simply as states. 

The number of electrons and the total spin can vary among the states. The time depen- 



ds! 



dence of the occupational probability Qa is given by the masters equation, [42, 

-^f- = X! O^P^aQp — ^a~*(iQa) , (5) 

where T^^p is the transition rate from state |a > to state \(3 >. 

The steady state solutions are obtained from the masters equation using = 0. The 
masters equations then mutually relate the occupational probabilities of various states. In 
addition, the occupational probabilities are normalized, that is, the sum of all occupational 
probabilities is equal to one. The masters equation in the steady state and the normalization 
condition are sufficient to determine the occupational probabilities. The current through the 
left barrier is obtained as 

/L = |e|EEr^-/3Q«, (6) 

a f3 

where T^^^ is the contribution of the left lead to the transition rate Fq^/j, taken with a 
positive or negative sign depending on weather the transition gives a positive or negative 



contribution to the current, respectively. 
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43| 



Because the tunnel density of states in the ferromagnets is spin-dependent, Qa will depend 
on the relative magnetic orientations of the leads. We set the magnetization of the left 
ferromagnetic lead to be always up. The magnetization of the right ferromagnetic lead can 
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be up or down, corresponding to values of parameter cr: if cr = 1 the magnetizations are 
parallel and if a = — 1 the magnetizations are antiparallel. 

The tunnel-densities of states at the Fermi level of spin up and spin down electrons in the 
left lead are N-^ = Nav{l + P) and A^j^ = Nav{l — P), respectively, where A^^^ is the average 
tunnel-density of states at the Fermi level per spin band, and P is the spin-polarization 
of the ferromagnet. Similarly, the tunnel densities of states at the Fermi level in the right 
lead are N-^ = Nav{l + crP) and Ni = Nav{^ — crP) for spin- up and spin-down electrons, 
respectively. 

We assume that electron spin is conserved in the tunnel process across a single tunnel 
junction. This assumption is justified by the fact that samples without nanoparticles have 
a large and weakly voltage dependent TMR. |6;] In that case, the spin-up (down) bands 
in the leads are tunnel-coupled only to the discrete levels in the nanoparticle with spin up 
(down). This is also valid because the single-electron states in the nanoparticle have well 
defined spin, since the spin-orbit coupling in the nanoparticle is weak, as discussed earher. 

The tunnel rate between the leads and a discrete level i is proportional to the tunnel- 
density of states in the leads and can be written as rj(l -|- P) for spin- up electrons and 
rj(l — P) for spin-down electrons, where Fj is referred to here as the bare tunnel rate. 

Next we use the masters equation to obtain the current in the regime defined in Sec. . 
To summarize, at zero bias voltage the number of electrons on the nanoparticle is even. We 
investigate the region of the IV-curve where only one extra electron can be added to the 
nanoparticle, within the first step of the Coulomb-staircase. This is the bias-voltage region 
where the saturation of the spin-polarized current with bias voltage is observed. We consider 
highly asymmetric junctions in resistance, Rl ^ Rr- A positive bias voltage is applied on 
lead R relative to lead L, so first an electron tunnels into the nanoparticle from the left lead, 
across the high-resistance junction, and after that the nanoparticle discharges one electron 
into the right lead via the low resistance junction. Finally, the internal spin-conserving 
relaxation rates are much larger than both tunnel rates, z^(co') ^ F-^ ^ F^. 

In this regime, it is convenient to divide the states into four groups. Group 1 contains 
the states which do not have an added electron and which are fully relaxed with respect 
to spin-conserving transitions. We describe these states in detail before discussing other 
groups. 

Consider the states which do not have an added extra electron displayed in Fig. [3l The 
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internal relaxation of the states in the figure must involve a spin-fiip process. So the states 
are fully relaxed with respect to spin-conserving transitions, thus they are in group 1. In 
the figure, the charging energy has been added to the single-electron energy levels of the 
nanoparticle, to make it clear which unoccupied single-electron energy levels are accessible 
for tunneling in from the left lead. As seen in the figure, the states from group 1 can be 
labelled \i >, where i = —N, —N + 1, N — 1, N, where = 4 for the example shown in 
the figure. 

is the number of unoccupied single-electron states above the Fermi level of the nanopar- 
ticle into which an electron can tunnel in from the left lead. It is related to the bias voltage. 
If the level spacing 6 is constant, then A^ = e c^^Cr ~^cb) / where Cl and Cr are the 
capacitances of the left and the right tunnel junction, respectively, and Vcb is the Coulomb 
blockade threshold voltage. {4Q] 

A^, the number of single-electron states that are accessible for tunneling-in, is different 
from the number of single-electron states that can discharge an electron (M). For example, 
if the nanoparticle is in state |0 >, then an electron from the left lead can tunnel only into 
one of the unoccupied single-electron levels labelled 1,2,...,A^ with either spin. After an 
electron tunnels in, the nanoparticle can discharge either the added electron or any other 
spin-up or spin-down electron from one of the doubly occupied single-electron states labelled 
0, -1, -2, -M (M = 12 in the figure), consistent with Coulomb blockade. 43| If the 
level spacing is constant, then M = eVcs/S. 

Group 2 contains states which do not have an added electron and which can internally 
relax by spin-conserving transitions. Group 3 contains states with an added extra electron 
and which are fully relaxed with respect to spin-conserving transitions. Finally, group 4 
contains the states with an added extra electron and which can internally relax by spin- 
conserving transitions. 

The next task is to determine the steady state occupational probabilities. Using //(cu) ^ 
YR ^ max(r'^, u^^luj)), several approximations can be made in the masters equation. 

First, it is shown in the appendix that the occupational probabilities Qi of the states from 
the first group are larger than the occupational probabilities of the states from the second, 
third, and the forth group, by a factor of and respectively. The occupational 

probabilities of the states from groups 2-4 can thus be neglected in the normalization con- 
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FIG. 3: Many-electron states from group 1. Black-filled and white-filled circles indicate occupied 
and unoccupied single-electron states, respectively. Only one electron can tunnel into the nanopar- 
ticle at a time because of the Coulomb Blockade; the displayed states are found before an electron 
tunnels into the nanoparticle from the left lead and after an electron tunnelled out of the nanopar- 
ticle into the right lead. The states are fully relaxed with respect to internal spin-conserving 
transitions. 



dition, so 



N 



i=-N 



with an error of order ^ -C 1, in the limit ^{uj) 3> r^(cc;). 

Second, if 1/(0;) 3> F^ 3> max(F^, //•^-^(o;)), the occupational probabilities of the states 
from the second, third, and the forth group of states can be eliminated from the masters 
equations in an explicit way. The elimination leads to a set of linear equations that relate 
the occupational probabilities within the space of states Gi. These equations are refereed 
to here as the renormalized masters equation: 



\ ^ -pren 



-pren 



(7) 



where T\^o is the renormalized transfer rate from state |m > to state |i >. 

The transfer \m >^ \i > takes place either directly, via a spin-flip transition, or indirectly 
via intermediate states. In the leading order of 7^,f^; ^m^i is equal to F^^j (direct 
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transition rate from m to i) plus the sum of the rates of transitions from state |m > into 
the intermediate states, weighted by the probabihty of transfer from the intermediate states 
into the state |i >: 

rzZi = r„^i + E rm^an(a ^ i), (8) 

a 

Il{a — > i) is the probabihty that the nanoparticle in intermediate state jo; > will transfer 
into state \i >. 

In the following, we will obtain F^ff^ in an intuitive way. This approach emphasizes 
understanding and enables one to obtain the occupational probabilities Qi from the renor- 
malized masters equation directly, without explicitly solving the masters equation. In the 
appendix, we will derive r^!!^j from the masters equation and show that the intuitive ap- 
proach is accurate within controlled approximations. 

NANOPARTICLE SPIN PROBABILITY DISTRIBUTION IN ABSENCE OF SPIN- 
RELAXATION 

In this section we examine the regime where no spin-relaxation takes place, ^ ^ 
^ v^^{uj). In that case, all transfers between the states within the space of states \i > 
are indirect. They involve a sequential electron tunneling process, in which the nanoparticle 
goes into an intermediate state with an extra added electron. r^!!^j is equal to the sum of 
the rates of transitions from state |m > into the various intermediate states with an added 
extra electron, weighted by the probability of transfer from these intermediate states into 
the state \i >. 

Consider the nanoparticle in state \i >, which has spin ih, {i > 1). After an electron 
tunnels into the nanoparticle, the nanoparticle is in an intermediate state with spin {i+l/2)h 
or (i — l/2)h. Then, after an electron tunnels out, the nanoparticle spin changes again by 
h/2 or —h/2, so the final spin after a sequential tunneling process can be {i — l)h, ih, and 
(i + 1)^. So the renormalized transfer rate from state \i > into states |m > is nonzero only 
ifm = i — lorm = iorm = i + l. 

The nanoparticle will transfer from state |i > to state \i + 1 > indirectly, if a spin-up 
electron tunnels in from the left lead and then a spin-down electron tunnels out into the 
right lead. A spin-up electron can tunnel into any of the unoccupied single-electron states 
k with spin-up, k — i + l,i + 2, ...N. After tunneling in, the nanoparticle instantly relaxes 

16 



via spin-conserving transition. A spin-down electron can then tunnel out from any of the 
occupied single-electron states j with spin down, j = —M,—M + 1, Then, if the 
nanoparticle after tunneling out is left in an excited state, it will relax instantly into state 
\i + l>. 

The rate for this spin-up-tunnel-in / spin-down-tunnel-out sequential process is obtained 
from Eq. [HI by summing over all the intermediate states with an added extra electron, 

N 

r^:^+l= E r,(i + p)7r,(i) (9) 

k=i+l 

where 7rj(|) is the probability that the nanoparticle with an added spin-up electron will 
discharge a spin-down electron. 

The rate at which a spin-down electron discharges is J2J1-m ^fi^ ~^P)i where is the 
bare tunnel rate between level j and the right lead. Similarly, the rate at which a spin-up 
electron discharges is Y^^^_m Tf{l + aP). The probabilities and vrj(J,) are proportional 
to the spin-up and spin-down discharge rates, respectively. The total discharge probability 
is one, T^i{]) + T^i{[) = 1, so we obtain 

^^(i) = :..r..Jl^ (10) 



1 + 



^—i] = — M 3 

As discussed previously, in our model the number of discharging levels M ^ 1. We make 
a further assumption that M ^ N, which is valid not too far from the conduction threshold 
for sequential electron tunneling through the nanoparticle. In this case, the denominator in 
Eq. [10] is close to one, within a factor of i/M ~ N/M, and we obtain 

l-aP 1 + aP 



TT. 



.a) = ^^,vr,(T) 



2 2 

The approximation N/M <C 1 enhances the probability to discharge a spin-down electron 
at the expense of suppressing the probability to discharge a spin-up electron. So the ap- 
proximation increases the spin accumulation efficiency. This approximation is not essential 
for our model to work, but it simplifies further calculations. 

Substituting into Eq. [9l we obtain 

rir^+i = A+liP), (11) 

where 

^,P).(i±£)ii^f:r.. (12) 
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Now we consider the indirect transfer \i >^ >• In this transfer process, an electron 
with spin-down tunnels in from the left lead, followed by a discharge of an electron with 
spin-up into the right lead. Following an analysis similar to the above, we find Tl^^_^ = 
A_,+i(-P). The left hand side of Eq. [7] becomes Qi{A+i{P) + A_i+i(-P)). Using Eqs. [H] 
and [23 the steady state Eq. [7] becomes 

Qi(A,+i(P) + A_,+i(-P)) 

= Qi-iAiP) + g,+iA_,(-p), (13) 

where i = 1,2, ...N. At i = N, in Eq. [13] we must put y4iv+i(P) = and Qn+i = 0. 

A similar analysis to the above leads to a steady state equation for the states with negative 
nanoparticle spin \ — i >, i = 1,2, N: 



g_,(A,+i(-p) + A_,+i(p) 

= Q-^+iA{-P) + Q_,_iA_,(P), (14) 

where i = 1,2, ...N. Ati = N, in Eq. [H we must put An+i{-P) = and Q-n-i = 0. The 
final equation necessary for finding Qi is Y.iL-NQi = 1- 

We calculate the probability distribution numerically for arbitrary A^, Fj, and Qi, allowing 
for fluctuations among different tunnel rates Fj across the left tunnel junction. Once the 
probability distribution is determined, the current though the nanoparticle is calculated 
from Eq. ^ which can be shown to be: 



T N J N 

FI i=-N 1^1 j=l 

N / N i \ 

T.qA E 2r,+ E r,(i-p) + 

i=l \j=i+l j=-i+l / 

N / N i \ 

E 2r,+ E F,(i + p) (15) 

i=l \j=i+l i=-i+l / 

To illustrate these equations, we plot the nanoparticle states along an axes, as shown in 
figure [H We indicate various transition rates obtained from Eqs. [TT] and [261 For simplicity 
we assume that all bare tunnel-in rates are the same (Fj = F). In that case 

rr-,.^(iv-„(l±^i^r, (16) 
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FIG. 4: Rates of tunnel and internal relaxation transitions between nanoparticle states |z > in 
Fig. [3l for the antiparallel magnetic orientations. For simplicity, the bare tunnel rates to the left 
lead (Ffc = F) are assumed to be independent of single-electron energy levels k. F+ = F(l + P)^/2, 
and F_ = F(l - Pf/2. In the parallel magnetization configuration, substitute r+ and F_ with 
F' = F(l — P^)/2. The wavy line indicates spin-relaxation processes. 



and 



rr_:,_, = (N + ,)^2^f)%^r. (17) 



The renormalized master equation is 

[{N - i)T{l + P)(l - aP)/2 + {N + t)T{l - P)(l + aP)/2] 

= Qi-i{N - 1 + i)r(i + p)(i - <jP)/2 + 

Qi+i ((iV + z + l)r(l - P)(l + (tP)/2) , (18) 
and the current through the nanoparticle, from Eq. [T5|, becomes 

/ = |e| {2NV - 2 < i > PT) , (19) 

where < i >= J^^^.^iQi- 

One notices from Fig. H] that the renormalized rate of transition \i >— |i + 1 > decreases 
as — i, when i increases from zero to N. Similarly, the renormalized rate of transition 
\i >—>■ \i — 1 > increases as N + i, when i increases from zero to N. So, if i is near N, the 
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rate of|i>— — l>is much larger than the rate of \i >^ \i + 1 > and in the steady 
state, Qi must be a rapidly decreasing function of i when i is near A^. Similarly, Q^i must 
be rapidly decreasing with i when i approaches A^. 

In the parallel magnetization orientation (cr = 1) the transition rates are symmetric 
around i = 0, as explained in Fig. HI In that case Qi has a maximum at i = 0, indicating 
that there is no spin-accumulation in the nanoparticle. This is a situation similar to spin- 
accumulation in large systems. 

In the antiparallel configuration of the leads, the rate of transition \i >— > ji + l > is 
proportional to (1 + P)^ and the the rate of transition \i >— \i — 1 > is proportional 
to (1 - P)2. In that case, transition rate \i |i+l>is larger than transition rate 

|i + 1 >^ \i > when i is close to zero, leading to spin-accumulation in the steady state. 
The spin accumulation is limited because the rate of transition \i \i + 1 > decreases in 
proportion with N — i and the rate of transition \i >— 1^ — 1 > increases in proportion with 
N + i, as discussed above. 

In the following two sections we focus on limit ^ 1, where we find the analytical 
solution of Qi and /. In the analytic solution we assume that the tunnel rates Fj across the 
left tunnel junction are independent of i, Fj = F. 



Spin-polarized current at large bias, N ^ 1 

At the maximum of Qi (the mode), which is found at 2 = io, the following condition is 
satisfied: F[-,^,^ = F^",^,^. This leads to 

2PN , , 

^0 = (20) 

in leading order of A^ in the antiparallel magnetic configuration [p = —1). Thus, in the most 
probable state of the nanoparticle, the chemical potential of spin-up electrons is shifted up 
by 2^0(5 relative to the chemical potential of spin-down electrons, 

= TTi^- (21) 

Next we obtain the fluctuations around the mode. We make a conjecture that the prob- 
ability distribution Qi is sharply peaked around the maximum ai i = io, so that the 
fluctuations around io are weak compared to A^. Our numerical calculations show that 
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rms{i) ~ vA^, confirming the conjecture. We can write i = io + j and expand Qi around 
i = Iq-. Qi = Q + Q'j + Q"j'^/2 + ... . Substituting into Eq. [TS], we obtain a differential 
equation, in leading order of V^, 

This linear differential equation can be solved analytically. Using a boundary condition 
Q{-^N) ^ 1, the normalized solution is 

/ 



1 + 

Qi ~ 1= exp 



2 ' 



(23) 



which is a Gaussian distribution with fluctuation 

/ATI -p2 



rms[i] 



2 1 + P2 • 

A similar analysis leads to the probability distribution in the parallel magnetic configu- 
ration: 

^-T^^^^W)- (24) 

and /|| = |e|2A^r, from Eq. [H 

Fig. [5] displays the electron distribution function in the nanoparticle in the parallel and 
antiparallel magnetization configurations. The difference in chemical potentials of spin-up 
and spin-down electrons in the antiparallel state is proportional to A^, the number of energy 
levels available for tunneling-in, according to Eq. [2T1 

Spin accumulation in the nanoparticle is well defined if the relative fluctuation, 

1 _p2 

rms{i)/io 



2PV2N 

is smaller than 1. For example, if P = 0.1, must be > 12 in order to have a well defined 
spin accumulation. If < 12, the time dependence of the nanoparticle spin will exhibit 
significant noise. In large systems where the level spacing is negligibly small, A^ is typically 
^ 1 and the fluctuations are thus negligible. 

Using Eq. [19] and the results of this section, the tunnel magnetoresistance can be shown 
to be 

P < ?■ > 2P^ 

™^ - T - T^- (2«) 
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FIG. 5: A and B: Probability distribution of electrons in parallel and antiparallel magnetization 
configuration, when the number of levels and M are large. 



4J] The Julliere's value is larger than the TMR predicted 



which is the JuIIieres's formula, 
theoretically. But this theory is valid in a different regime than ours and our approxima- 
tion slightly enhances the spin-accumulation efficiency, as discussed earlier. Another recent 
numerical calculation [? ] in a regime similar to our own, and an assumption of infinite spin 
lifetimes confirms the JuUiere value as a maximum of the TMR in a nanoparticle. 

NANOPARTICLE SPIN PROBABILITY DISTRIBUTION WITH SPIN- 
RELAXATION 

If spin- relaxation rate is not negligible, then the nanoparticle in state \i > can undergo 
a spin-flip transition into state \i — 1 >. The nanoparticle spin changes by ^ in a spin-flip 
transition, hence there are no other final states in the space of states Gi after a spin-flip 
transition. The renormalized transfer rate is 



prert 



A_,+i(-P)+(], 



(26) 



where fli is the total transfer rate from |z > to |2 — 1 > that includes an internal spin flip 
transition. This transfer can take place directly or via intermediate states. Starting from a 
state \i > in Fig. [3l an electron occupying the highest single-electron energy level i with spin 
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up can make a transition into the lowest unoccupied single-electron energy level + 1 with 
spin down. In this process the energy difference is {2i — 1)6 and the rate is i'^^[[2i — 1)6]. 
This is a direct transition. 

Alternatively, starting from the same sate \i >, an electron occupying the highest single- 
electron energy level i with spin-up can make a transition into the next-to-the lowest un- 
occupied single-electron energy level —i + 2 with spin-down. The nanoparticle is left in 
an excited state, which is an intermediate state. This is followed by a spin-conserving re- 
laxation transition, which is instantaneous in our model, and the nanoparticle ends in the 
state 1^ — 1 >• Similarly, starting from state \i >, an electron occupying the next-to-the 
highest single-electron energy level i — 1 with spin-up can make a transition into the low- 
est unoccupied single-electron energy level —i + 1 with spin-down, leaving a hole. This is 
followed again by an instantaneous spin-conserving relaxation transition, which brings the 
nanoparticle into the state |i — 1 >. The energy differences for the spin-flip processes are 
the same, {2i — 2)6, assuming equal level spacing. The total spin-flip transition rate with 
this energy difference is 2i'sF[{2i — 2)6]. 

Taking into account all spin-flip processes with varying energy differences, we find the 
spin-relaxation rate 

a = E 0' + ^>sf[{2i - j - 1)6]. (27) 

3=0 

In general, we expect to be rapidly increasing with u: u^^ {oj) ~ u^. In that case 

VLi increases with energy faster than v^^iu): Vli ~ i""*"^. 

In this analysis we neglect higher order spin- flip transitions \i >— \i — 2 >, \i >— >■ 
|i— 3 >,... , where the environment would receive angular momentum 2h, 3h,... , respectively. 
For the low energy states of the nanoparticle, we expect that the probability of the higher 
order processes to be much smaller than the probability of the first order processes. 

If we assume that spin-relaxation process is governed by phonon emission and Elliot- 
Yafet scaling, then Qi = i'^^{6) = ave-ph{6). In that case, the spin-relaxation rate increases 
very rapidly with the excitation energy. From Eq. [271 we find VL2 = 46fii, ^3 = 371f2i, 
fli = 1596^1, and 

a ~ 1.6aue-phi6)i^ = l.en^i^. (28) 
In the following we will show that the rapid increase of Qi with i causes the saturation 
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behavior of the spin-polarized current. It is coincidental that the spin-relaxation rate Qi 
of the nanoparticle increases with fifth-power of the excitation energy, analogous to its 
temperature dependence in bulk. 

Spin-relaxation in the nanoparticle reduces the spin-accumulation. Consider again fig- 
ure m The difference between rates and r_ moves the distribution mode from zero to 
positive i, as discussed in the previous section. The spin relaxation rate Qi is added to F^, 
reducing the difference between the upward and downward rates. So the distribution mode 
moves downward relative to the mode at f2j = 0. 

The distribution mode io is obtained from V^^^li^ig = Tl^Hi^i^. It follows that in the 
antiparallel state, in leading order of A^, the mode satisfies the following equation 



We set the spin-relaxation rate f2i to be much smaller than the tunneling rate F, fii ^ F, 
because if Qi ^ F, the spin-polarized current would be close to zero, contrary to our 
measurements. In that case = 

N[2P/{1 + P2)] for small N. 
As N increases, both iq and the spin-relaxation term increase. The spin-relaxation rate 
Qi increases with i much more rapidly than i (as i""'"^). The spin- relaxation begins to reduce 
the spin-accumulation when the two terms become comparable, around 



As A^ increases further, the spin-relaxation term becomes dominant, and the mode is ob- 
tained from fijo = A^[2P/(1 + P^)]. Assuming fij ~ i"+2 follows that at large A^, which 
corresponds to large bias voltage, ~ A^i/("-+2)_ xhis is a much weaker dependence on A^ 
than N[2P/ (1 + P^)]. So there is a crossover in spin-accumulation versus N, from linear de- 
pendence to a much weaker dependence. Similarly, the spin-accumulation crosses over from 
liner V-dependence into a much weaker V-dependence at large V, because N, the number 
of single-electron energy levels available for tunneling in, is linear with bias voltage V, as 
discussed before. 

One consequence of the rapid increase of Qi with i is an asymmetric spin probability 
distribution. Our numerical calculations show that the probability that the nanoparticle 
spin is below the mode is larger than the probability that the nanoparticle spin is above 
the mode. In that case, the average nanoparticle spin is smaller than the most probable 




(29) 




N 



P 
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nanoparticle spin, < i >< io. Nevertheless, as long as the width of the distribution is much 
smaller than the mode, one can substitute io for < z > in Eq. [19] and the spin polarized 
current I-^^ — = 2\e\P < i > T ^ 2\e\PiQ exhibits the same crossover with bias voltage. 
This explains the saturation of spin-polarized current with bias voltage observed in our 
experiments. 

The crossover condition can be rewritten as = NPT. So, at the crossover point, 
the rate of the spin-flip transition with energy difference Afi (Eq. [2T|) corresponds to the 
spin-polarized current, e.g. Q{Afi) ^ (/|| — I^)/\e\. 

Now we discuss the spin-relaxation effects assuming that the relaxation is mediated by 
phonon emission, in accordance with the Elliot- Yafet relation Eq. [2S1 Figure EJ-A displays 
the distribution function at large bias voltage, = 50, obtained by numerical calculations 
for P = 0.1. Qi/T is varied from 0.25 to 0.25 ■ 10"''. Also shown is the Gaussian distribution 
in Eq. [231 which is valid in absence of any spin-relaxation. At Qi = 0.25-10"^ the probability 
distribution is very close to the Gaussian, indicating that the spin relaxation is negligible. 
As Qi/T increases, the mode shifts downwards and the distribution becomes asymmetric. 
At Qi/T = 0.25, the mode is located at io = 1, showing that there the spin accumulation is 
very weak. 

The distribution mode is now obtained from 

2P 1.6Qiif. , , 

'^TTpJ = '» + r(TTi%- 

At large N, when the spin relaxation term dominates, 

/2PA^r\^/^ 

that is, the mode crosses over from being linear with to being proportional to N^l^ . 

Fig. [6]-B displays spin-polarized current versus A^ obtained by numerical calculations, 
for different fii/F. The crossover from linear to a much weaker dependence is evident for 
fii/F = 0.25 and 0.25 ■ 10~^. The crossovers in io and — /|| are equivalent, as discussed 
before, so at large A^, — ~ N^l^ . 

At fii/F = 0.25, the rate of spin-relaxation with energy difference b is small compared to 
the tunneling rate. But r22/F = 11.5, so the spin-relaxation rate with an energy difference 
> 2b is large compared to the tunneling rate. As a result, the spin-polarized current crosses 
over already at the second single electron energy level. In particular, we find that the 
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FIG. 6: A: Probability distribution functions of the nanoparticle many-electron states in the an- 
tiparallel state, when the number of energy levels available to tunnel in is = 50. B: Spin-polarized 
current A/ = — versus N, the number of energy levels available to tunnel in. and bias 
voltage are related linearly. 

contribution from the second single-electron energy level to the spin polarized current is 
approximately a 3rd of the contribution from the first single-electron energy level. 



FITTING 



To illustrate how the model derived in this paper can be used to interpret measurements 
of the spin-polarized current, we discuss sample 1 in Ref. 6[. In Fig. [7]-A we display the 
I-V curve at T = 0.03K, obtained in the parallel magnetization configuration and in the 
regime where the electron discharge rate is much faster than the tunnel in rate, as required 
by our model (we have reversed the sign of bias voltage compared to that in Ref. 6|). The 
I-V curve increases in discrete steps at voltages corresponding to discrete energy levels _of 
the nanoparticle. The average tunnel in rate is obtained as F =< 61 > /2\e\ = 1.5MHz, 
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FIG. 7: A: I-V curve in sample 1 near the sequential tunneling threshold. B: Circles: Spin polarized 
current A/ = /|| — I = I^i versus A^. Line is the best fit to the model. 



'Ti 



where < 5/ > is the average current step. The average level spacing S = O.SmeV corresponds 
to a spherical nanoparticle with diameter 6nm. The electron g- factors are very close to 2, 
confirming weak spin-orbit scattering. 

The spin polarized current, AJ = — I^, versus N is shown by circles in Fig. [3-B. At 
small A^, energy levels are well resolved, and A^ is observed directly as shown in the figure. 
At large bias voltage, where energy levels are broadened, we find A^ as ec{V — Vcb)/^, 6| 
where c is the capacitance ratio that converts from voltage to nanoparticle energy and Vcb 
is the Coulomb blockade voltage threshold. 

We fit AJ obtained from the model versus A^, using a fixed tunnel rate F = 1.5MHz and 
two fit parameters, P and Qi. We assume Qi ~ l.Gi^Qi, which is valid when spin relaxation 
is mediated by phonons following Elliot- Yafet scaling, as discussed earlier. The best fit 
parameters are P = 0.23 and Qi = I.IF ^ 1.6MHz. This value of the spin-relaxation rate 
is consistent with what we estimated based on qualitative discussions in Ref. (gI 

The best fit does not fully saturate with bias voltage, as seen in Fig. [6]-B and [Tj-B, but 
exhibits a crossover to a weaker dependence with A^. Nevertheless, we find the agreement 
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between the model and the data to be quahtatively good. An outlying point at = 3 in 
Fig. [3-B is attributed to the large tunnel rate via the third energy levels as seen in Fig. O-A, 
because we use a constant tunnel in rate in the fit. Such a large tunnel in rate arises from 
the natural statistical fluctuations of the tunnel rates among different levels. 

If we assume that the spin- relaxation rate is energy independent, Qi = Qi, the spin- 
polarized current in our model will have a linear dependence with A^. In that case, not 
only does the best fit have five times larger chi-square, but also the best fit parameters 
are un-physical, P = 0.045 and Qi = 0.067r. In particular, the measured spin-polarized 
current via the low energy levels would be 25 times larger than the maximum theoretical 
value corresponding to TMR = 2P^/(1 + P^) at P = 0.045. Consequently, an energy 
independent spin-relaxation rate cannot explain our results, and the energy dependence of 
the spin-relaxation rate plays a critical role in interpreting the spin-polarized current through 
the nanoparticle in the regime of well resolved energy levels and weak spin-orbit scattering. 

Our model neglects higher order spin relaxation processes. We suggest that there could 
also be transitions where two electrons flip their spin simultaneously, emitting a phonon with 
angular momentum 2h. The probability of such a transition path must be much weaker 
relative to the probability of the 1st order process (weaker by a factor of a) . But with 
increasing i, the number of transition paths behind a second-order process would increase 
much faster with i than in Qi, leading to a larger exponent in the dependence of the spin- 
relaxation rate on i. As a result, at high excitation energies, second order or even higher 
order spin-relaxation process can dominate spin-relaxation, leading to another crossover in 
A/ versus V, weakening the V dependence further. Thus, inclusion of the the higher order 
spin-relaxation processes at high excitation energies would improve the agreement between 
model and measurements. 

CONCLUSION 

In conclusion, we have derived a simple model for calculating spin-accumulation and 
spin-polarized current via discrete energy levels of a metallic nanoparticle in the regime of 
weak spin-orbit scattering. The model is valid in a relatively narrow range of sample param- 
eters. However, a large percentage of samples fabricated by lithography have parameters 
in that range. The energy dependence of the spin-relaxation rate causes a signiflcant sup- 
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pression of the bias voltage dependence of the spin-polarized current at large bias voltage. 
In particular, if the spin-relaxation rate increases with excitation energy as a power of n, 
u^^luj) ~ cu", then the spin polarized current at large bias voltage V will increase propor- 
tionally to \/i/{«+2)^ which is a dependence much weaker than linear. The crossover between 
the linear dependence at low voltage and the much weaker dependence at large voltage oc- 
curs when the the spin-polarized current/ |e| is equal to the spin- relaxation rate with the 
energy difference given by the spin-accumulation. The model leads to the spin-relaxation 
rate 1.6MHz in an Aluminum nanoparticle of diameter 6nm, for a transition with an energy 
difference of one level spacing. 

This work was performed in part at the Georgia- Tech electron microscopy facility. This 
research is supported by the DOE grant DE-FG02-06ER46281 and David and Lucile Packard 
Foundation grant 2000-13874. 

APPENDIX: MASTERS EQUATION 

In the steady state the occupational probabilities are time-independent and the masters 
equation [5] becomes 

Qa 5Z ^ce^P = X! Ql3^l3^oi- (31) 

The occupational probabilities satisfy the normalization condition Yl,aQa = 1- 

In this section we reduce the problem of solving Eq. [3l] into a simpler problem. This is 
done by identifying a group of states which have much larger probability than the states 
outside the group and then eliminating the occupational probabilities of the states outside 
the group. 

We use inequalities Va^p ^ max(r^, i^a^/s). Va^p is the spin-conserving relaxation 

rate from state \a > to state |/5 > given by EqUJ in which u = E^—Ep is the energy difference 
between the initial and the final state. We calculate at zero temperature, so Ua-tp = if 
Ea < Ej3. Va^i3 is uouzcro only if the initial and the final state have the same spin. In 
this regime the electron transport through the nanoparticle is in equilibrium with respect 
to spin-conserving relaxation. 

J-'a^p is the spin-flip relaxation rate between the states. It can be calculated using Eq. [2l 
^a^B is nonzero only if Sg — = ±h. 



29 



In the analysis, we separate the states that are fully relaxed with respect to spin- 
conserving relaxation processes from those that are not. The latter states decay at rate 



the number of electrons is even (before tunneling in after tunneling out) or odd (after tun- 
neling in before tunneling out). From Eq. [HI] it follows that the occupational probabilities 
of states that are not fully relaxed with respect to spin-conserving transitions are strongly 
suppressed in the limit Ua-^p ^ ^ max(r^, v^^p), because the occupational probability 
of those states have u^^p in the denominator. 

We are considering only the states with one extra added electron, as discussed in sections 
and . Let Qi denote the occupational probabilities of the states without an added extra 
electron that are fully relaxed by spin-conserving transitions. These states are from group 
Gi and are shown in Fig. [31 Let Pj denote the occupational probabilities of the states with 
an added extra electron that are fully relaxed by spin-conserving transitions; these states 
are from group 6*3. 

Next we consider the excited states. Q'^ denotes the occupational probabilities of the 
states without an added electron that are not fully relaxed by spin-conserving transitions; 
these are the states from group G2- Finally, P'^ denotes the occupational probabilities of the 
states with an added extra electron that are not fully relaxed by spin-conserving transitions; 
these are states from group 6*4. 

The masters equations are next written for the states in each group. For groups 1 and 2, 
the equations are 



Va^Pi whereas the former states decay at a rate max (F^, v^^) or 



F , depending on weather 




(32) 



Q'a I]mGGi(^( 



a 



■m 



a- 
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+ Ej'GGs Ta^j' + E/3'gG4 



F. 
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+ 2j'6G3 Pj'^j'^a + J2f3'£G4 Pp'^P'^a, (33) 

where i E Gi, a E G2. Rates Fj^j/, Fj^^/, Ta^j/, and Fq^^/ on the left hand sides are 
proportional to the tunnel-in rate F'^ across the left junction. Similarly, rates Fj/^j, F^/^j, 
Ta^jf, and F^/^^ on the right hand sides are proportional to the tunnel-out rate F^ across 
the left junction. 

For groups 3 and 4, the equations are 



Pi 



E 



(34) 



P' 



+ Ej'GG3('^/3-j' + '^l^j') + E/3'GG4('^/3-/3' + 



5F > 
/3^/3'J 



^j'€G3 ^i'^i'^p + E/3'eG4 Pp'i^P'^P + ^p^-^p)^ 



(35) 



where j G G3, /3 G G4. 

First we estimate the orders of magnitude of the occupational probabilities. The tran- 
sition rates are assumed roughly constant to within an order of magnitude and pulled out 
of the sum, and Q\ P, and P' refer to the probability of finding the particle in a 
state from group 1, 2, 3, and 4, respectively. 47|] Entering the orders of magnitude of var- 
ious terms in Eq. [331 we obtain Q'{iy + u^^ + F^) ~ Qu^^ + Q'{u + z/^^) + PT^ + P'F^, 



where v and v^^ indicate orders of magnitude of Va^jj and 



SF 
a^p-i 



respectively. In the limit 



z/ > F^ > max(F^, z/^^), one gets Q'v ~ Qv^^ + (P + P')T^. 

Similarly, substituting orders of magnitude in Eq. [32], one obtains Qiy^^ + V^) ~ Qv^^ + 
Q'{v + u^^) + PF^ + P'F^. In the limit z/ > F^ > max(F^, z/"^^), one gets Q{u^^ + F^) ~ 
Q'v + (P + P')T^. Combining the two order of magnitude estimates we obtain 



,SF 



+ F^ 



^rnM^^^^ « g. 



(36) 
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Hence, the occupational probabilities Q'^ of the excited states from group G2 are much 
smaller than the occupational probabilities Qi of states that are fully relaxed with respect 
to spin-conserving transitions, an expected result. 

Next, we substitute into Eq. [33 the orders of magnitude of various terms and obtain 
P'(r^ + z/ + z/^^) QT^ + Q'T^ + Pv^^ + P'{v + v^^). Using z/>r^ >max(r^,z/^^) and 
Eq. I36[ one obtains P'u ~ QV^ + Pv^^. Finally, Eq. [M] leads to order of magnitude estimate 
P(r^ + z/ + v^^) ~ QV^ + QT^ + Pv^^ + P'{v + z/^^), which leads to PV^ ~ QT^ + PV. 
Combining the estimates, we obtain 

p ~ « g,p' ~ p^J^ « p (37) 

So, the occupational probabilities of states with an added electron and fully relaxed with 
respect to spin-conserving transitions are much smaller than the occupational probabili- 
ties of states without an added electron and fully relaxed with respect to spin-conserving 
transitions, as expected from the large asymmetry in tunnel resistance. In addition, the 
occupational probabilities of excited states with an added electron are much smaller than 
the occupational probabilities of states with an added electron and fully relaxed with respect 
to spin-conserving transitions. 

Using the estimates in Eq. [36] and [33, the leading order of the masters equation for 
groups 1 and 2 are 



SF 



a' (iG2 i^a 



(38) 



Q'a EmGGi ^a^m + Ha'£G2 ^ol^ol' 
- SmeGi Qml^m-^a a' eG2 Q'a'^a'-,a 



(39) 



where i E Gi, a E G2. For groups 3 and 4, the equations are 
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^31 



(40) 



Pp T>j'€G3 ^P^i' + S/3'eG4 '^/S-/?' 

= LmeGi + I]/3'eG4 Pf3'^l3'^l3, (41) 

where j & G^, f3 & G4. The relative errors of the terms in these equations are smaller than 
the maximum of max(r^, z/'^'^)/i/, Tl/Tji, u^^/Tji, and T^/v. 

Now we begin the process of elimination. The first step is to eliminate the excited 
states. In principle, the space of excited states is very large, if the excited states include 



multiple electron hole pairs. 



33] These multiply excited states are generated by tunneling if 



the particle remains excited longer than the time of a sequential tunneling cycle, so in our 
limit the occupational probabilities of the excited states with multiple electron-hole pairs 
are negligibly small. In this regime, the space of excited states G2 can be restricted to the 
excited states that can undergo a direct spin-conserving transition into a state from Gi. 
Similarly, the space of excited states G4, can be restricted to those excited states that can 
undergo a direct spin-conserving transition into a state from G3. 

Consider Eq. HH which represents an equilibrium condition for a state |/3 >, /5 G G4. For 
a given state \j >, j G Gs, we select a subspace within G4, such that the nanoparticle in a 
state from the subspace can relax via spin-conserving transition into the state \j >. That 
is, for any state \P > within the subspace, 1^/3^ j 7^ 0, and for any state |/3 > outside the 
subspace, i^/3^j = 0. Then we sum Eq. [41] over the subspace, which leads to 

J2'l3eG4 ^j'&G-i Pp^P^r + E^GG4 S/3'GG4 P'p^P^P' = 
I^mGGi I^/3eG4 Q-rrXm-*(5 + I]/3gG4 '^(i'&Gi Pp'^P'^Pi 

where I]/3gG4 sums only over the states within the subspace, so that vp^j 7^ 0. 

In the second term on the left hand side of this equation, where (3 is restricted within the 
subspace, /?' becomes also restricted within the subspace, because of the spin-conservation 
(since states \[3 > and \j > have the same spin, and states \P' > and \P > also have the same 
spin, must be restricted within the subspace). The second term becomes Yl'p^p'eGi Pp^fS^fS'- 
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Similarly, in the second term on the right hand side, fp'-rp is nonzero only if states \ j3' > and 
1/5 > have the same spin as the spin in state \j >, and the term becomes J^'is^p'eGi Pp'^/s'-'is- 
Exchanging the indices f3 and the second terms on the left hand side and the right hand 
side cancel. 

Now consider the first term on the left hand side. As stated above, only those states \(3 > 
that can undergo a spin- conserving transition into the state \j > are in the sum over (3. It 
follows that z//3_^j/ is nonzero only if j' = j, since varying states \j' > from space G3 have 
different spin. The first term on the left hand side becomes J2'i3eG4, ^'p^P-^j- ^^"^ remove 
the prime and sum instead over the entire space G4, because automatically restricts 
the sum to the subspace. One obtains 

The left hand side in this equation is the same as the second term on the right hand side 
of Eq. Ho] and thus it can be eliminated, which leads to 

Next we perform a similar elimination process using Eqs. [38] and [39] For a given i G Gi, 
we sum Eq. [39] over the excited states |q; >, a G G2, for which Ua-^i 7^ 0, 

Z]aGG2 SmGGi Q'a^a-*m + S0GG2 Sa'GG2 Qa^a^a' = 
J2aGG2 SmGGi Qm^^m^a + SjGGs J2aGG2 Pj^j^'^'^ 

Following a similar analysis to that above, the second term on the left hand side is equal to 
the third term on the right hand side, and the first term is nonzero only for m = i, so one 
obtains 

Substituting this equation and Eq. [12] into Eq. [38], and following several lines of alge- 
bra, one obtains the renormalized masters equation in the space of states Gi, with the 
renormalized rate 



pren 

m— >j m— >i •<-^aGG2 m— >a 
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X .. ^-^^--^a- r (43) 

This expression is the same as that obtained intuitively in the main text. The right hand 
side in the first row is the renormahzed spin-flip rate, ^m^i = I'^^i + Y^'a&cf^ ^m^a- ^rn-^i 
is the direct spin-flip transition rate and the sum is taken over the excited states that can 
undergo a spin-conserving transition into the flnal state \i >. 

In the spin-flip process the spin decreases by Ti, hence i^^^i 7^ only if AS'^ = S** — S*^ = 
-h, and z/^l„ ^ only if = 5^ - 5f = -fi. It follows that ^ only if z = m - 1. 

In that case, one can denote fim-^m~i = = i^m^m-i + EagGT"'^" ^m^a- in this 
equation is the renormahzed spin-flip rate given by Eq. [271 

Next we examine the second and the third row of Eq. |l3l The sum over j is taken over 
intermediate states with an extra added electron. Writing the rates T^^y explicitly (which 
is not shown here), it can be seen that the rates of tunneling in to the various unoccupied 
single electron states are multiplied by the rates of tunneling-out from the various occupied 
single electron states. The contribution to r^!!,j is nonzero only if z = m — 1 or i = m or 
i = m + 1, because a sequential tunneling process changes the spin by —1 or or 1. 

Consider the expression rj^!!^^_^]^ in Eq. HHl In the second row term J2jeG:i(Xm^j + 
Z]^'gG4^° Tm-^/s) X {■■■), j indicates the state obtained by adding an electron into the lowest 
unoccupied single-electron level with spin up. The sum over /3 is taken over the excited 
states that can relax via spin-conserving transition into \j >. One obtains J2jeG3{^m~>j + 
Y^'pecf'^ '^m-^p) ^ (■••) — J2k=m+i^k{^ + -P) X (...), whcrc Ffc is the tunnel-in rate across the 
left lead into an unoccupied single-electron energy level k. This expression is the same as 
that in Eq. E 

The third row in Eq. |13] can be interpreted as the probability that the nanoparticle in 
state j will discharge a spin-down electron, which is the same as 7rj(|) in Eq. [91 Substituting 
the tunnel rates into Eq. [431 explicitly, one obtains that the third row of Eq. [l3l is the same 
as Eq. M 

In summary, the model of electron transport from the intuitive approach is derived in 
this appendix, using the masters equations. 
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